!***************************************
finish
/clear
/filn, APDL_Example  
/title, vortex
/units,si                 !ùʵλ

!*****************************************
pi=3.14159265
rr=8.75e-3                !߻Բľ
rb=4.2176e-3              !߻Բ뾶
alaf=30                   !߷
hw=0.03                   !ݸ
wo=0.0835                 !̰뾶
h=0.02                    !̺
exx=2.07e11               !òϵĵģ
nuxyqq=0.3                !òϵĲɱ
gxyqq=0.5*exx/(1+nuxyqq)  !òϵļеģ
densqq=7850               !òϵܶ
alpqq=9.5e-6              !òϵϵ
layer1=6                  !Ԫ

!****************ǰ***************************
/prep7  
!************趨Ԫͼ***************
et,1,95                   !1ŵԪ-solid95
mp,ex,1,exx               !趨ϵĵģ
mp,gxy,1,gxyqq            !趨ϵļеģ
mp,nuxy,1,nuxyqq          !趨ϵĲɱ
mp,dens,1,densqq          !趨ϵܶ
mp,alpx,1,alpqq           !趨xϵ
mp,alpy,1,alpqq           !趨yϵ
mp,alpz,1,alpqq           !趨zϵ


*afun,deg                 !ǺöǵĴС

!ͨؼ
!ȡȦڲཥߵĹؼ
*DIM,inner,ARRAY,111,2                                 !111x2Ķά
*do,ii,0,1100,10                                       !趨ѭii0仯1100Ϊ10
inner(ii/10+1,1)=rb*(cos(ii)+(ii+alaf)*pi/180*sin(ii)) !ȦڲཥߵĹؼx
inner(ii/10+1,2)=rb*(sin(ii)-(ii+alaf)*pi/180*cos(ii)) !ȦڲཥߵĹؼy
*enddo                                                 !ѭ
!*status,inner,1,111,1,2,1                             !бʾ

!ȡȦཥߵĹؼ
*DIM,outer,ARRAY,108,2                                 !108x2Ķά
!ע⣬ȦıںȦڡཥʼƫƫΪ30ȣii30ʼ仯
*do,ii,30,1100,10                                      !趨ѭii30仯1100Ϊ10
outer(ii/10-2,1)=rb*(cos(ii)+(ii-alaf)*pi/180*sin(ii)) !ȦཥߵĹؼx
outer(ii/10-2,2)=rb*(sin(ii)-(ii-alaf)*pi/180*cos(ii)) !ȦཥߵĹؼy
*enddo                                                 !ѭ


!**********************ģ************************
!****************֪Ŷؼ********************
*do,ii,0,1100,10                                       !趨ѭii0仯1100Ϊ10
k,ii/10+1,inner(ii/10+1,1),inner(ii/10+1,2)            !ڲཥ߹ؼ
*if,ii,ge,30,then                                      !жϵѭii=30֮
k,ii/10+201,outer(ii/10-2,1),outer(ii/10-2,2)          !ཥ߹ؼ
*endif                                                 !жϽ
*enddo                                                 !ѭ


!************⿪߲䶨3㣬Աߵķ***********
*do,fai,30,50,10/3
x2=rb*(cos(fai)+(fai-alaf)*pi/180*sin(fai))         !߷õx
y2=rb*(sin(fai)-(fai-alaf)*pi/180*cos(fai))         !߷õy
k,fai*3/10+192,x2,y2                                !ؼ
*enddo                                              

!*************ɵ㹹*********************************
*do,aak,1,106,5                                     !趨ѭΪ5ʼΪ1ֹ106
bsplin,aak,aak+1,aak+2,aak+3,aak+4,aak+5            !ڲཥ
*enddo                                              !ѭ

*do,aal,201,306,5                                   !趨ѭΪ5ʼΪ201ֹ306
bsplin,aal,aal+1,aal+2,aal+3,aal+4,aal+5            !ཥ
*enddo                                              !ѭ

!************ߣΪ׼*******************
*do,jj,1,111,5                                      !趨ѭΪ5ʼΪ1ֹ111
l,jj,jj+200                                         !ߣΪ׼
*enddo                                              !ѭ



!************̽ߵĩɵ**************
lfillt,22,67,0.001
lfillt,67,44,0.001
!*************ɵһ߶*******
lcomb,68,67                       !67ԭ߶45Ϊһ߶                     
lcomb,69,67                       !46ԭ߶45Ϊһ߶   
l2tan,-1,23                       !߶233еֱ
ldel,45

!**********************߹****************
*do,abc,2,22,1 
al,abc,abc+45,abc+22,abc+44
*enddo
al,23,68,1,46

!**********ɾĹؼ**********************
allsel                      !ȫѡ
lsla                        !ѡϵ
ksll                        !ѡϵĵ
ksel,invert                 !ѡ
kdele,all                   !ɾ

allsel
*get,maxnum,kp,,num,max     !ؼ

!*************ɵz****************
k,maxnum+1                    !ؼ
k,maxnum+2,,,hw               !ؼ㣨ϵxyһ£zΪhw
l,maxnum+1,maxnum+2           !ͨ㶨߶
lesize,45,,,layer1,2          !¶߶趨ַ 
vdrag,all,,,,,,45             !¶ֱ

NUMMRG,all                    !Ա
NUMCMP,all                    !ѹ


cylind,0,wo,,-h,,360          !
allsel                        !ȫѡ
vsel,u,volu,,23               !ȥ
vglue,all                     !е

!*******************************
type,1                        !趨Ԫ1
mat,1                         !ָ1
esize,0.006,0,                !趨ֵԪС
vsel,all                      !ѡе
vsweep,all                    !ӳ񻮷

!!!!!!!!!!!!!ڵ֮佨Լ*******
vsel,s,volu,,23               !ѡе
eslv                          !ѡиڵ̵ĵԪ
nsel,u,ext                    !unselectڵ̵ԪϵĽڵ
nsel,r,loc,z,0                !λZ0reselectڵ
ceintf,,all                   !ѡԼ
/PBC,ALL,,1                   !ʾʩӵԼ
finish